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We study the entanglement preservation of two qubits locally interacting with their reservoirs. We 
show that the existence of a bound state of the qubit and its reservoir and the non-Markovian effect 
are two essential ingredients and their interplay plays a crucial role in preserving the entanglement 
in the steady state. When the non-Markovian effect is neglected, the entanglement sudden death 
(ESD) is reproduced. On the other hand, when the non-Markovian is significantly strong but the 
bound state is absent, the phenomenon of the ESD and its revival is recovered. Our formulation 
presents a unified picture about the entanglement preservation and provides a clear clue on how to 
preserve the entanglement in quantum information processing. 

PACS numbers: 03.65.Yz, 03.67.Mn 



I. INTRODUCTION 

Entanglement is not only of fundamental interest to 
quantum mechanics but also of great importance to quan- 
tum information processing [1]. However, due to the in- 
evitable interaction of qubits with their environments, 
entanglement always experiences degradation. Entangle- 
ment sudden death (ESD), a phenomenon in which the 
entanglement between two qubits may completely disap- 
pear in a finite time, has been predicted theoretically [2[ 
and subsequently been verified experimentally [U, indi- 
cating specific behavior of entanglement that differs from 
that of coherence. From the point of view of applications, 
ESD is apparently disadvantageous to quantum informa- 
tion processing. 

Recently, Bellomo et al. [4( found that the entangle- 
ment can revive after some time interval of ESD and thus 
extends significantly the entangled time of the qubits. 
This remarkable phenomenon, which has been experi- 
mentally observed is physically due to the dynam- 
ical back action (that is, the non-Markovian effect) of 
the memory environments [HQ- However, in many cases 
the finite extension of the entangled time is not enough 
and thus it is desired to preserve a significant fraction of 
the entanglement in the longtime limit. Indeed, it was 
shown Q that some noticeable fraction of entanglement 
can be obtained by engineering structured environment 
such as photonic band-gap materials [1, Q . According to 
these works, it is still unclear whether the residual entan- 
glement is fundamentally due to the specific structured 
materials or due to certain physical mechanisms. Is there 
any essential relationship between ESD and/or its revival 
phenomena and the residual entanglement? 

In this work we focus on these questions and elucidate 
the physical nature of the residual entanglement. Before 



* anjhong@lzu.edu.cn 
t phyohch@nus.edu. sg 



proceeding, it is helpful to recall the physics of quan- 
tum electrodynamics of a single two-level atom placed 
in a dielectric with a photonic band gap [1, The 
coupling between the excited atom and electromagnetic 
vacuum in the dielectric leads to a novel photon-atom 
bound state Q, in which the fractional atomic popula- 
tion on the excited state occurs, also known as population 
trapping (Toj . This result has been verified experimen- 
tally for quantum dots embedded in a photonic band-gap 
environment fill ]. The population trapping has been di- 
rectly connected to the entanglement trapping due to the 
structured environment [?|. Here we reveal that there 
are two essential conditions needed to preserve the en- 
tanglement. One is the existence of the bound state be- 
tween the system and its environment, which provides 
the ability to preserve the entanglement, and the other is 
the non-Markovian effect, which provides a way to pre- 
serve the entanglement. Our result can reproduce ESD 
2] when the non-Markovian effect is neglected. The phe- 
nomenon of ESD and its revival discussed in Ref. [3] 
results from the non-Markovian effect when the bound 
state is not available. The interplay between the avail- 
ability of the bound state and the non-Markovian effect 
can lead to a significant fraction of the entanglement pre- 
served in the steady state. We verify these results by 
considering two reservoirs modeled by the supcr-Ohmic 
and Lorentzian spectra, respectively. The result provides 
a general method on how to protect the entanglement by 
engineering the environment. 

This paper is organized as follows. In Sec. [Til the 
model of two independent qubits in two reservoirs is in- 
troduced. By exploring the eigen-spectrum of the model, 
we derive the condition for the formation of bound states 
between the qubits and their respective reservoirs and 
discuss its profound consequence on the dynamics of the 
open two-qubit system. In Sec. IIII1 the entanglement 
preservation caused by the formation of a bound state 
and the non-Markovian effect is studied numerically in 
different cases of the environmental spectral density. Fi- 
nally, we end with a short conclusion and a discussion 
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about the experimental feasibility of our result in Sec. 
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II. MODEL AND DECOHERENCE DYNAMICS 

We consider a system consisting of two independent 
subsystems, each of which contains a qubit coupled to a 
zero-temperature reservoir. Due to the dynamical inde- 
pendence between the two subsystems, we can first in- 
vestigate the single subsystem, then extend our studies 
to the double-system case. The Hamiltonian of each sub- 
system is [12| 

H = uj a + a- + u k b\ b k + (g k cr+ b k + g% <r_ b[ ) , (1) 
fe fe 

where a± and loq are the inversion operators and transi- 
tion frequency of the qubit, b\ and bk are the creation and 
annihilation operators of the fcth mode with frequency ui k 
of the reservoir, and g k denotes the coupling strength be- 
tween the atom and the radiation held. 

To check the spectrum of the Hamiltonian we hrst solve 
the eigenvalue equation 

H\<p E )=E\<p E ). (2) 

If only one excitation is present in the system at 
zero temperature initially, then \lpe) — co|+,{Ofc}) + 
Sfc c fe l — j Ifc): where |±) is the atomic excited (or ground) 
state, and |{0fc}) and \l k ) are the vacuum state and the 
state with only one photon in the fcth mode of the reser- 
voir, respectively. Substituting Eq. ((T|) into Eq. @, one 
has 

V(E) ee loq f°° ^-du> = E, (3) 
Jo w - E 

where J(ui) = Ylk \dk\ 2 S(u; — oj k ) is the spectral density 
of the reservoir. The solution of Eq. ([3]) highly depends 
on the explicit form of J[uS). If the reservoir contains 
only one mode to', then J(ui) — g 2 5(w — ui'). This is the 
ideal Jaynes-Cummings (JC) model 12], in which two 
bound states in one excitation sector are formed and as 
a result the dynamics of the system displays a lossless 
oscillation. When the reservoir contains infinite modes, 
one can model J(w) by some typical spectrum functions 
such as the super-Ohmic or Lorentzian forms. 

We hrst consider the super-Ohmic spectrum J(lu) = 

3 — / 

?7^e _w/Wc , where 77 is a dimensionless coupling constant 

and lu c characterizes the frequency regime in which the 
power law is valid 13]. It corresponds to the case in 
which the reservoir consists of a vacuum radiation field, 
where g k oc ^fUJk [l2| . The existence of a bound state re- 
quires that Eq. ((3]) has at least a real solution for E < 0. 
It is easy to check that the solution always exists if the 

condition y(0) < (i.e. ojq — 2??— | < 0) is satisfied. 
Otherwise, no bound state exists. This condition can 



be fulfilled easily by engineering the environment. For 
the Lorentzian spectrum it is found that a criterion for 
the existence of a bound state cannot be obtained an- 
alytically. In this case one can use the diagrammatic 
technique shown later. 

The existence of a bound state has a profound implica- 
tion on the dynamics of the single-qubit system such as 
the inhibition of spontaneous emission [8ll9l.ll4|. Further- 
more, it also has an important impact on the entangle- 
ment dynamics of a two-qubit system, which is governed 
by the master equation [la ] 
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P(t) = ^{-ifl(t)K<r!,pi+r(t)[2a^(t)a? 

n=l 

-<<r»p(i) - p(t)a$a»]}, (4) 

where O(t) = -Im[^g], and T(t) = -Regg]. It is 
shown that co(t) satisfies 

C (t) + i0J Q C (t) + f Co(T)f(t-T)dT = 0, (5) 

Jo 

where f(t - r) = J °° J "(w)e~ ltJ ( t - T ) duj . The time- 
dependent parameters f2(t) and T(t) play the role of 
Lamb-shifted frequency and decay rate of the qubits, re- 
spectively. The integro-differential equation ([5J contains 
the memory effect of the reservoir registered in the time- 
nonlocal kernel function f(t — r) and thus the dynamics 
of the qubit displays a non-Markovian effect. If f(t—r) is 
replaced by the time-local function, then Eq. (|4]) recov- 
ers the conventional Born-Markovian master equation, 
where the parameters become constants [l6[, that is, 

r - 7tJ(w ), fio = wo - P/ °° 4^? with P denoting 
the Cauchy principal value. 

The dynamical consequence of the bound state is de- 
coherence suppression [10|. If a bound state is absent, 
then Eq. ([3]) has only complex solutions. Physically this 
means that the corresponding eigenstate experiences de- 
cay from the imaginary part of the eigenvalue during the 
time evolution, which causes the excited-state popula- 
tion to approach zero asymptotically and the decoher- 
ence of the reduced qubit system. While if a bound state 
is formed, then the population of the atomic excited state 
in the bound state is constant in time because a bound 
state is actually a stationary state with a vanishing de- 
cay rate during the time evolution. So there will be some 
residual excited-state population in the long-time limit. 
Due to this decoherence suppression, we expect that the 
formation of a bound state plays a constructive role in 
entanglement preservation under the non-Markovian dy- 
namics, as shown in the following section. 

III. ENTANGLEMENT PRESERVATION 

To study the entanglement dynamics of the bipar- 
tite system, we use the concurrence to quantify entan- 
glement [T7|. The concurrence is defined as C(p) — 
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FIG. 1. (Color online) Entanglement dynamics of the two 
qubit system with local super-Ohmic reservoirs, (a) Diagram- 
matic solutions of Eq. ((3]) with different parameters. C(t) as 
a function of time is shown in (b): (u) c ,rf) = (Q.Tujo, 0.2), (c): 
(oj c ,v) = (0.7oj , 1.0) and (d): (<jJ c ,v) = (3.0u;o, 0.2). The pa- 
rameter q is taken as 0.7. For comparison, C(t) under the 
Markovian approximation has also been presented by using 
the same parameters. 



max{0, y/X\ — VA2 — VA3 — V^-TK where the decreasing- 
order-arranged quantities Xi are the eigenvalues of the 



matrix p(<jy ® (jy)p*{<jy <g> ). Here p* means the com- 
plex conjugation of p, and <r y is the Pauli matrix. It can 
be proved that the concurrence varies from for a separa- 
ble state to 1 for a maximally entangled state. Consider 

an initially entangled state |-0(O)) = a\ ) +/3|++), 

where \a\ 2 + |/3| 2 = 1. Then C(t) can be calculated as 
C(t) = max{0,(5(t)}, where 

Q(t) = 2|«/3||c (£)| 2 - 2|/3| 2 |c W| 2 [l - |c (t)| 2 ], (6) 

which indicates that the time-dependant factor of the 

excited state population (|co(i)| 2 ) determines solely the C'(t) = max{0,2e 
entanglement dynamics. 

Now we are ready to study the entanglement dynam- 
ics of the two-qubit system. Consider first the super- 
Ohmic case. Figure [1] shows the entanglement dy- 
namics in different parameter regimes [i.e., {uj Cl ri) — 
(0.7w ,0.2),(0.7cj ,1.0) and (3.0w , 0.2)]. For the first 
two parameter sets the bound state is absent, while 
for the last one it is available, as shown in Fig. QJa). 
Whether the bound state exists or not plays a key role 
in the entanglement preservation in the longtime limit. 
When the bound state is absent, the residual entangle- 
ment approaches zero in a long enough time, as shown by 
the non-Markovian lines in Figs. (Hb) and (He). The dif- 
ference between these two cases is that Fig. Hfb) shows 
the weak-coupling regime, where the non-Markovian ef- 
fect is weak, while Fig. QJc) shows the strong-coupling 
regime, where the strong non-Markovian effect leads to 
an obvious oscillation. When the bound state is avail- 
able, the situation is quite different, as shown in Fig. 



FIG. 2. (Color online) The decay rate V(t) as a function of 
time in the non-Markovian and Markovian cases. The param- 
eters used are uu c — S.Oujq an d r\ — 0.2. 



QJd). The entanglement first experiences some oscilla- 
tions due to the energy and/or information exchanging 
back and forth between the qubit and its memory en- 
vironment [lH |, then approaches a definite value in the 
longtime limit, where the decay rate approaches zero af- 
ter some oscillations, as shown in Fig. [21 The entan- 
glement preservation is a result of the interplay between 
the existence of the bound state (providing the ability to 
preserve the entanglement) and the non-Markovian ef- 
fect (providing a way to preserve the entanglement). The 
claim can be further verified by the fact that the entangle- 
ment preservation is absent in the Markovian dynamics, 
as shown by the Markovian lines in Figs. Hfb), [Tfc) and 
QJd), where the entanglement displays sudden death ir- 
respective of the availability of the bound state. This is 
because the Markovian environment has no memory and 
the energy and information flowing from the qubit to its 
environment is irreversible and the decay rate remains 
fixed (see Fig. In this case one has 
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which shows a finite disentanglement time when \a\ < 
|/3 1 [2j. In short, different to the need for a structured 
environment as emphasized in Ref. Q, our discussion 
clearly reveals two essential conditions to preserve the 
entanglement: the availability of the bound state and 
the non-Markovian effect. 

This discussion focused on an almost maximally entan- 
gled initial state by taking a = 0.7. In Fig. Owe show 
the results for different initial states with different initial 
entanglement. With decreasing initial entanglement, the 
residual entanglement also decreases in the longtime limit 
and, finally, ESD happens for a — 0.3. The result can be 
understood from Eq. ©. On the one hand, the residual 
entanglement is determined by co(oo), which is directly 
related to the property of the bound state. On the other 
hand, the residual entanglement is also determined by 
the competition between the first and the second terms 
in Eq. ([B]), which is dependent of the initial state. 
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FIG. 3. (Color Online) The residual entanglement for dif- 
ferent initial states with a — 0.7,0.5 and 0.3. The other 
parameters used are the same as those in Fig. [T] 



In order to make a comparative study and confirm our 
observations we consider the Lorentzian spectrum if the 
reservoir is composed of a lossy cavity, 
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(8) 



where 7 is the coupling constant and A is the spectrum 
width. This model has also been studied in Ref. Q, 
where the lower limit of the frequency integral in f(t — r) 
was extended from zero to negative infinity. This exten- 
sion is mathematically convenient but the availability of 
the bound state is missed. Here we follow the original 
definition of the frequency integral ranges. 
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FIG. 4. (Color online) The entanglement dynamics with the 
Lorentzian spectrum, (a) Diagrammatic solutions of Eq. Q 
with different parameters A = O.lwo, 2.0o;o and I5cjo. (b) C(t) 
as a function of time for the corresponding three parameter 
regimes. The insert in (b) shows the decay rate as a function 
of time. The other parameters used are 7 = 3.0cjo and a — 
0.7. 

Our model with the Lorentzian spectral density corre- 
sponds exactly to the extended damping J-C model flij . 
It is noted that the strong coupling of the J-C model 
has been achieved in circuit QED [l9( and quantum dot 
[2fj| systems. Figure [4] shows the entanglement dynamics 
of the qubits under the Lorentizan reservoir for differ- 
ent spectral widths in the strong-coupling regime. When 
A = O.Iwq, Eq. ([3]) lacks the bound state. According 



to our discussion, there is no residual entanglement in 
the longtime limit. This is indeed true, as shown in Fig. 
|4jb). However, it is noted that before becoming zero 
the entanglement exhibits "sudden death" and revives 
for several times. This is an analog of the central result 
found in Ref. |4j, that is, the phenomenon of ESD and 
revival. Apparently, this is due to the non-Markovian ef- 
fect with the revival being a result of back action of the 
memory reservoir. The situation changes with increasing 
the spectral width, and the bound states become avail- 
able. A significant fraction of the entanglement initially 
present is preserved in the longtime limit, where the de- 
cay rates shown in the insert of Fig. |4jb) approach zero 
in these cases. Likewise, the physical nature of the en- 
tanglement preservation is still the interplay between the 
bound state and the non-Markovian effect. The stronger 
the coupling is, the more striking the entanglement os- 
cillates as a function of time and, consequently, the more 
noticeable the non-Markovian effect is, as shown in Fig. 
[5] For 7 = 0.2wo, the system is in the weak-coupling 
regime, where the bound state is also not available. As 
a result, ESD is reproduced in this case. 
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FIG. 5. (Color online) The same as Fig. [4] but A = 15.0wo is 
fixed and 7 = 0.2oio, 2.0wo, and S.Oujq. 

In Fig. [5] we present a phase diagram of the entan- 
glement in the steady state for the Lorentzian spectral 
density. In the large-7 and small-A regime, the system 
approaches the J-C model. In this situation the strong 
back action effect of the reservoirs makes it difficult for 
the qubit system to form a steady state. The entangle- 
ment oscillates with time but has no dissipation. In the 
small-7 and large-A regime, the non-Markovian effect is 
extremely weak and our results reduce to the Markovian 
case. In the limit of a flat spectral density, the Born- 
Markovian approximation is applicable and the system 
has no bound state. This is the case of ESD Q. 



IV. CONCLUSIONS AND DISCUSSION 

In summary, we have studied the entanglement protec- 
tion of two qubits in two uncorrelated reservoirs. Two es- 
sential conditions for preserving the entanglement are ex- 
plored: the existence of the bound state of the system and 
its reservoir and the non-Markovian effect. The bound 
state provides the ability of the entanglement preserva- 
tion and the non-Markovian effect provides the way to 
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FIG. 6. The phase diagram of the residual entanglement in 
the steady state for the Lorentzian spectrum. 



protect the entanglement. Previous results on the entan- 
glement dynamics in the literature can be considered as 
specific cases where these two conditions have not been 
fulfilled at the same time. The result here provides a 
unified picture for the entanglement dynamics and gives 
a clear way on how to protect the entanglement. This is 
quite significant in quantum information processing. 

The presence of such entanglement preservation ac- 
tually gives us an active way to suppress decoherence. 
This could be achieved by modifying the properties of 
the reservoir to form a bound state and to approach the 
non-Markovian regime via the potential usage of reservoir 



engineering [2l| - |23l |. Many experimental platforms (e.g., 
mesoscopic ion traps cold atom BECs [22j, and the 
photonic crystal materials [3| have exhibited the con- 
trollability of decoherence behavior of relevant quantum 
systems through optimally designing the size (i.e., modi- 
fying the spectrum) of the reservoir and/or the coupling 
strength between the system and the reservoir. It is also 
worth mentioning that a proposal aimed at simulating 
the spin-Boson model, which is relevant to the one consid- 
ered in this paper, has been reported in a trapped ion sys- 
tem [24{ ■ On the other hand many practical systems can 
now be en gin eered to show the novel non-Markovian ef- 
fect [!, [25M28j . All these achievements demonstrate that 
the recent advances have paved the way to experimen- 
tally simulate the paradigmatic models of open quantum 
system, which is one par t of the newly emerging field of 
quantum simulators [29j. Our work sheds new light on 
the way to indirectly control and manipulate the dynam- 
ics of a quantum system in these experimental platforms. 
It provides a clue to preserving the entanglement in quan- 
tum information processing. 
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